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The heavy-fermion system CeCue-^Au^ exhibits a quantum critical point at Xc ~ 0.1 separating 
nonmagnetic and magnetically ordered ground states. The pronounced non-Fermi-liquid behavior at 
Xc calls for a search for the relevant quantum critical fluctuations. Systematic measurements of the 
inelastic neutron scattering cross section S{q, ui) for a; = 0.1 reveal rod-like features in the reciprocal 
ac plane translating to two-dimensional (2d) fluctuations in real space. We find 3d magnetic ordering 
peaks for x = 0.2 and 0.3 located on these rods which hence can be viewed as 2d precursors of the 
3d order. 
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Continuous quantum phase transitions which occur in 
a strict sense only at temperature T = are driven by 
quantum fluctuations instead of thermal fluctuations as 
for ordinary classical phase transitions This leads 

to unusual and rich behavior even at finite temperatures 
in the neighborhood of the critical point. Because of the 
uncertainty principle the energy scale of fluctuations in- 
troduces a time scale which leads to an intricate coupling 
of static and dynamic critical behavior. For instance, 
the critical behavior of the specific heat will depend on 
the dynamical critical exponent z relating the typical 
lifetime and correlation length ^ of critical fluctua- 
tions, ^T- oc Such a quantum phase transition can be 
achieved by changing a coupling parameter which plays 
a role analogous to temperature in ordinary phase tran- 
sitions. In recent years, many physical realizations of 
quantum phase transitions have been found. The case 
of a magnetic-nonmagnetic transition in heavy-fermion 
metals is particularly interesting because of the involve- 
ment of itinerant electrons. 

Excitations of a system of interacting itinerant elec- 
trons in a metal, i.e. quasiparticles, are usually de- 
scribed within the Fermi-liquid theory, with the spe- 
cific heat C (X T, a Pauli susceptibility independent of 
T, and an electrical resistivity contribution Ap oc 
due to quasiparticle-quasiparticle scattering. Interac- 
tions renormalize the quasiparticle masses with respect 
to the free-electron mass toq- Even in heavy-fermion sys- 
tems with quasiparticle masses as high as several 100 
mo, Fermi- liquid behavior is the rule rather than the ex- 
ception 1^. In heavy-fermion systems, the coupling pa- 
rameter tuning the magnetic-nonmagnetic transition is 
the (antiferromagnetic) exchange interaction J between 
4/ or 5/ magnetic moments and conduction electrons 
1^. If it is strong, a local singlet state is formed via 
the Kondo effect around each 4/ or 5/ site, leading 
to a nonmagnetic ground state. On the other hand, 
a weak (but non-zero) exchange interaction leads to a 



Rudermann-Kittel-Kasuya-Yosida coupling between mo- 
ments and hence to magnetic order. In the exemplary 
system CeCug-ajAua; doping of CeCue with the larger 
Au atom leads - via lattice expansion - to a weakening of 
the Kondo effect and hence to long-range antiferromag- 
netic order for x > Xc ~ 0.1, with a linear increase of the 
Neel temperature T/v oc {x — XcY, i.e. ^ = 1 @]. 

At Xc where T/v vanishes, i.e. around the quantum 
critical point, pronounced deviations from Fermi- liquid 
behavior occur. This non-Fermi-liquid (NFL) behavior 
is seen, e.g. in the specific heat where C /T oc — ln(r/To) 
over nearly two decades and in the resistivity where 
Ap oc T. It is precisely this NFL behavior at the quan- 
tum critical point that stirred a lot of interest ||] since 
it cannot be explained in terms of a transition driven by 
three-dimensional (3d) fiuctuations, because for an anti- 
ferromagnet with d — i and z — 2, C /T oc 1 — By/T and 
Ap oc r3/2 would be expected ||,|. 

A step forward towards the solution of the NFL puzzle 
in CeCu6-a:Aua; was to realize that 2d critical fiuc- 
tuations coupled to quasiparticles with dynamics will 
indeed lead to 7 = C/T oc -ln(T/ro), Ap oc T and 
/i = 1 as experimentally observed. Elastic neutron scat- 
tering experiments at 0.07 K on CeCu5.8Auo.2 with T/v = 
0.25 K showed, in addition to peaks attributed to (short- 
range) antiferromagnetic order, broad maxima along the 
a* direction that were much sharper in the h* direction 
0. This latter feature was interpreted in terms of fer- 
romagnetic planes perpendicular to the a direction (or- 
thorhombic notation) and thus provided a possible sce- 
nario of the d = 2, z — 2 universality class [Q. Without 
any direct evidence it is certainly hard to believe that 2c? 
correlations are dominating an intrinsically 3d alloy, even 
if the thermodynamics strongly support such a picture. 
Therefore it is essential to investigate the quantum crit- 
ical fluctuations directly by inelastic neutron scattering. 
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FIG. 1. a) g-scan along [1 0] in CeCus.gAuo.i at 
T = 50 mK. b) q-scans {ho I) in the a*c*-plane in the c* 
direction for fixed feo at T = 70 mK (all scans with fixed 
kf — I.ISA"^ and hoj — O.lmeV). The scans are shifted by 
150 counts to each other. 

The experiments were carried out at the triple-axis 
spectronieter IN14 at the Institut Laue-Langevin, Greno- 
ble with a fixed final neutron energy Ef = 2.7 meV 
{kf = 1.15A~^), giving an energy resolution (FWHM) of 
0.07 meV. The CeCug-xAua; single crystals were grown 
with the Czochralski method in a W crucible. The spe- 
cific heat of the sample with x = 0.1 exhibits the NFL 
behavior C/T (x — ln(r/To) as measured down to 60 mK, 
in agreement with previous samples of the same Au con- 
centration 0]. 

Fig. 1 shows q scans of the dynamic structure factor 
S{q^ w) in the reciprocal ac plane in two perpendicular 
directions at very low T < 100 mK at an energy transfer 
hco = 0.1 meV. The (h 0) scan (Fig. la) reveals a broad 
double maximum at (0.8 0) and (1.2 0). This double 
maximum is only resolved at small fuv. For instance, for 
hiu = 0.25 meV only a single broad feature centered at (1 
0) is seen [|j . Hence it may be thought of as developing 
from the broad maximum observed at the same q for fku 
= 0.3 meV in CeCug [§. Upon entering the magnetically 
ordered state for the x — 0.2 alloy, the double-peak struc- 
ture appears as a (quasi-) elastic feature for x ~ 0.2 that 
represents short-range ordering evidenced by a width in 
q that is considerably larger than the q resolution . 



Fig. lb shows that for x = 0.1 there is a very rich 
structure of S{q,u}) in the a*c* plane, as derived from 
scans {ho I) along c* for fixed h = ho. The peak at (1.2 
0) splits when moving away from the a* axis. The solid 
lines present Lorentzian fits with a width of (0.24 ±0.02) 
A^^ for all scans shown. However, the main point is 
that the peak height remains roughly constant across the 
whole Brillouin zone (cf Fig. 3a) . The width along c* is 
comparable to the width of the (1.2 0) maximum along 
{h 0) (cf. Fig. la). This suggests a rod-like feature of 
the dynamical magnetic response. It does not, however, 
extend along the a* axis as previously assumed |^] but in 
an oblique direction. 
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FIG. 2. Position of the dynamic correlations {x — 0.1, 
hu = O.lmeV, T < 100 mK) and magnetic Bragg peaks 
{x = 0.2 — 1.0) in the reciprocal ac plane in CeCu6-a:Aua; 
(data for x = 0.3 and 1.0 from @, for x = 0.5 from 
The vertical and horizontal bars indicate the line width for 
X — 0.1. The four rods are related by the orthorhombic 
symmetry (Pnma). Here we ignore a small monoclinic dis- 
tortion (< 1%) at low temperatures. The open symbols for 
X — 0.2 represent the short-range ordering peaks The in- 

set shows a resolution limited magnetic Bragg peak for x — 0.2 
at T = 50 mK (elastic scan at (1.375 /))■ 

Fig. 2 shows the peak positions in the a*c* plane for 
X — 0.1 and huj = 0.1 meV derived from Fig. 1 and from 
further measurements. To corroborate the rod-like na- 
ture of S{q,hoj — O.lmeV) for x — 0.1, further scans 
across the peaks were performed along independent di- 
rections, one in the a*b* plane along b* and one in the 
a*c* plane perpendicular to the rod-like feature. They, 
too, reveal a width of comparable magnitude as can be 
seen from Fig. 3b and c. In order to interpret S{q, lo) 
of CeCu5.9Auo.i, we recall that a rod-like feature in q 
space is related to a 2d correlation between Ce atoms 
in real space. The shaded rods in Fig. 2 can therefore 
be identified with planes in real space. These planes ex- 
tend into the b direction and into a direction in the ac 
plane given by next-nearest neighboring Ce atoms. Thus, 
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the observed quasi 2d correlations strongly support the 
proposed scenario [Q of 2d spin fluctuations coupled to 
quasiparticles with 3d dynamics, although the quasi 2d 
correlations are not ferromagnetic as initially supposed. 

The 2d fluctuations apparently are the precursor of the 
3d magnetic ordering. Indeed, the Bragg points for sam- 
ples not too far from the magnetic instability, e.g. x = 0.2 
and 0.3, are located on the rods for x = 0.1. For x = 0.2 
in addition to the rather broad double maximum at q = 
(0.8 0) and (1.2 0) we find resolution-hmited peaks 
at (0.625 0.275) and at lattice-equivalent positions in 
reciprocal space. The inset in Fig. 2 displays a (1.375 1) 
scan of such a Bragg peak. The main frame of Fig. 2 
shows that its position is indeed on one of the rods as 
is the position of the short-range order peaks along a*. 
However, we have not observed a 3d precursor for a; = 0.1, 
i.e. enhanced scattering intensity around the Bragg peak 
for X = 0.2. This is an important point in favor of the 2d 
scenario. For x = 0.3, the Bragg position remains almost 
unchanged while no short-range order peaks on the a* 
axis were detected. For x — 0.5 a sudden reorientation 
of the magnetic ordering vector is observed, with incom- 
mensurate order along a* with t = (0.59 0) ||ll| which 
is then roughly constant up to a; 1 ||l^. The reorien- 
tation of r occurring between x = 0.3 and 0.5 deserves 
further study. 

Returning to the quantum-critical point at x — 0.1, 
we recall that in the d — 2, z ~ 2 scenario we expect the 
following generic form of the ui- and g-dependent sus- 
ceptibility describing magnetic fluctuations in the plane 
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The imaginary part of the susceptibility is directly 
proportional to the magnetic structure factor ^(q, uj) 
measured with inelastic neutron scattering, S'(q, to) = 
(1 + nB{uj))lm.x^^{i^)f{q±), where ns(w) is the Bose- 
function. The smooth function f{qi_) describes the weak 
q-dcpcndcnce perpendicular to the planes, i.e. along the 
rod-like structures shown in Fig. |2[ q|| is the momen- 
tum in the plane, i.e. perpendicular the rods, qa and loq 
are constants that vary only slightly with temperature 
T and depend only weakly on the momentum along the 
rods. We expect that the above equation is valid (up to 
logarithmic corrections) for small momenta and frequen- 
cies WjWoq^/'Zo ^ ^bTk, where Tk ~ 6K is the Kondo 
temperature in the system. We have neglected the small 
anisotropy within the planes. The effective correlation 
length is given by ^. It is expected to vary strongly with 
T, wo/{qoO'^ ~ AksT with A being a constant of order 
one varying only logarithmically with T ]l2|,pt. 

The q-scans were performed with an energy transfer of 
0.1 meV (= 1.2K-fcB) which is large compared to the tem- 
perature of 70 mK. Therefore we expect that the width 
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FIG. 3. g-scans along (a) and across (b, c) the rod-like 
feature in CeCus.gAuo.i for huj — 0.1 meV and kf — 1.15 A^^ 
(T = 50/70 mK). There is only a weak g-dependence along 
the rods {q±), while transverse scans (gy) show peaks with 
nearly the same line width. 

of the peaks shown in Fig. |^ and ^ does not measure 
the true correlation length ^ but determines the ratio 



q^/ujo ~ (0.1 ± 0.02)(meV Ay 



(2) 



The important question arises whether the observed 
magnetic fluctuations can be related to the NFL be- 
havior of the thermodynamic quantities at the quan- 
tum critical point, i.e. to the logarithmically diverg- 
ing specific-heat coefficient. Actually, the prefactor of 
the specific-heat coefficient per area of a plane is fully 
determined by the quantum critical theory 7^'' = 
{n/12){qQ/ujQ)\n{TQ/T) where n is the number of spin 
components. We use n — 1, as the magnetic anisotropy 
Q suggests an Ising system. Tq is an unknown temper- 
ature scale of the order of Tk ■ To calculate the specific 
heat per volume one has to know the distance L of the 
planes - note that especially in an incommensurate struc- 
ture such a distance is only an effective quantity. Then 
the molar specific heat is given by 
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Vm is the volume per mol of CeCus.gAuo.i, the factor 
2 takes into account that the correlations show up in 
two different directions (see Fig. 2). This value has to 
be compared to the measured [|j specific heat coefficient 
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7 = 0.6 ^Jj^-i ln(To/T) which is indeed of the same or- 
der of magnitude as our estimate (^). From the crystal 
structure one would expect L to be of the order of 4 
to lOA while from (|) we obtain L « 2 — SA which is 
somewhat too small. However, one has to take into ac- 
count the considerable theoretical uncertainty, e.g. aris- 
ing from the definition of L in an incommensurate system 
or the unknown effective number of spin-components in 
this anisotropic system. In addition, the momenta and 
frequencies used in our analysis are quite large. 

Therefore, we think that the semi-quantitative agree- 
ment of the width of the rods in q space compared to 
the specific heat gives strong support for the idea that 
the 2d fluctuations are responsible for the observed NFL 
behavior. 
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FIG. 4. Energy scans at q = (1.375 0.275) in 
CeCus.gAuo.i for kf = 1.15 A^^ for different temperatures; 
o 70 mK, • 500 mK and ■ 5 K. For comparison an energy scan 
&t q = (1.8 0) for T = 70mK is displayed (□). The inset 
shows energy scans at (1.2 0) in CeCue-aiAui: for a; = 
(T = 80mK) and x = 0.1 (T = 50mK), kf = 1.15 A'^ The 
background is subtracted for all scans. 

As a final point, we discuss the energy dependence of 
the critical modes. Fig. 4 shows energy scans at exactly 
the q value of the magnetic order in CeCu5.8Auo.2 (cf. in- 
set of Fig. 2) . The solid lines indicate a fit comprised of 
an elastic Gaussian and a quasi-elastic Lorentzian (convo- 
luted with the resolution) with full-width F ~ 0.62,0.14 
and 0.11 meV for the latter at T = 5,0.5 and 0.07 K, re- 
spectively. The shift towards lu = with decreasing T is 
clearly seen at this q as opposed to the rather broad lu 
scan a.t q = (1.8 0) with F = 0.48 meV for T = 0.07 K 
(deduced from polycrystalline measurements). Due to 
the finite energy resolution and the large elastic back- 
ground it is difficult to obtain reliable values for F at 
low temperatures. For instance, assuming a weakly uj- 
dependent background reduces F at T = 0.07 K by a 
factor of 2-3. Nevertheless, the observed decrease of F 
seems to be slower than F = tJo/('i'oC)^ ^ ^ predicted for 



the d = 2, z = 2 quantum critical theory for T <C Tk ~ 
6K. It is important to note that the onset of 3cZ correla- 
tions would have the opposite effect at this q value where 
Bd magnetic order is expected. We note that in measure- 
ments of S{q,uj) concentrated on q = (1.2 0) a similar 
decrease of F towards low T with apparent leveling-off is 
observed iQ. On the other hand, the thermodynamics, 
notably the specific heat, shows a logarithmic increase of 
7 down to at least 0.06 K without a signature of a cor- 
responding energy scale. Therefore this point certainly 
requires further studies, e.g. on the role of disorder in 
our incommensurate system. 

In conclusion, we have identified the critical two- 
dimensional fiuctuations leading to non-Fermi-liquid be- 
havior in CeCus.gAuo.i with a systematic study of 
quasielastic neutron scattering. From the observed dy- 
namic susceptibility a semi-quantitative agreement with 
the prefactor of the logarithmic increase of the specific- 
heat coefficient is found. 
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Forschungsgemeinschaft . 
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